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^ i Abstract: The nonlinear mechanical response of carbon nanotubes, referred to as their 

g 2 "buckling" behavior, is a major topic in the nanotube research community. Buckling means 

^ s a deformation process in which a large strain beyond a threshold causes an abrupt change in 

d ■ 

g h the strain enregy vs. deformation profile. Thus far, much effort has been devoted to anal- 

5 ysis of the buckling of nanotubes under various loading conditions: compression, bending, 

^ ^ torsion, and their certain combinations. Such extensive studies have been motivated by (i) 

O . 

O 7 the structural resilience of nanotubes against buckling and (ii) the substantial influence of 

8 buckling on their physical properties. In this contribution, I review the dramatic progress in 

^ 9 nanotube buckling research during the past few years. 

on: 

Keywords: nanocarbon material; nanomechanics; nonlinear deformation 

^ : 

1. Introduction: Appeal of nanocarbon materials 

Carbon is a rare substance that takes highly diverse morphology. When carbon atoms form a three- 

. 5^13 dimensional structure, their glittering beauty as diamonds is captivating. When aligned in a two-dimensional 

14 plane, they make up just black graphite and lose their sparkle. In addition to these "macro"-scopic car- 

15 bon materials, several "nano"-carbon materials have been discovered in the past few decades, opening 

16 up new horizons in material sciences. It all began with the Ceo molecule (fuUerene), whose existence 

17 was predicted by Osawa [1] in 1970 and was discovered by Kroto et al. [2] in 1985. Subsequent studies, 

18 including those on carbon nanotubes by lijima [3] in 1991^ and on graphene by Novoselov et al. [4] 



^See Appendix A for the delicate issue about who should be credited as the first person to discover carbon nanotubes. 
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19 in 2004, have had a tremendous impact and driven developments in science and engineering around the 

20 turn of the century.^ 

21 Among the three types of nanocarbon materials, carbon nanotubes are attracting the greatest attention 

22 in both industry and academia. Research on carbon nanotubes has brought out two characteristics not 

23 usually seen in other fields. First and foremost is the sheer breadth of the research, which encompasses 

24 physics, chemistry, materials science, electronics, and fife science. The second characteristic is that 

25 basic research and applied research are extremely close to each other. A succession of phenomena of 

26 interest to scientists has been discovered like a treasure chest, each leading to an innovative application 

27 or development. Nowadays, it is difficult even for professionals in the nanotube research community to 

28 understand the progress being made outside of their field of expertise. 

29 One of the reasons why carbon nanotubes offer huge potential is the fact that mechanical deforma- 

30 tion causes considerable changes in electronic, optical, magnetic, and chemical properties. Thus, many 

31 studies on new technologies to utilize the correlation between deformation and properties are underway 

32 in various fields. For example, studies of nanoscale devices based on the change in electrical conductiv- 

33 ity or optical response resulting from deformation are one of the most popular trends in nanotechnology. 

34 Another important reason for nanotube research diversity is the concomitance of structural resilience and 

35 small weight, making realizable ultrahigh- strength materials for utilization in super-high-rise buildings 

36 and large aerospace equipment. Furthermore, applications of these low-density substances for aircraft 

37 and automobile parts will raise fuel efficiency and save energy, as well as dramatically reduce exhaust 

38 gas emissions and environmental impact. 

39 With this background in mind, I shall review recent development in a selected area of nonlinear 

40 mechanical deformation, the "buckling" of carbon nanotubes.^ Section 2 provides a concise explanation 

41 of the terminology of buckling, followed by a survey of different approaches used in nanotube buckling 

42 investigations. Section 3 details the two most interesting features observed during nanotube buckling 

43 process, i.e., the structural resilience and sensitivity of nanotube properties against buckling. Sections 4 

44 to 8 are the main part of this paper, illustrating nanotube buckling under axial compression (§4), radial 

45 compression (§5), bending (§6, 7), and torsion (§8). Section 9 presents a universal scaling law that 

46 describes different buckling modes of nanotubes in a unified manner. The article is closed by Section 

47 10 that describes several challenging problems whose solutions may trigger innovation in the nanotube 

48 research community. The list of references (over 190) is fairly extensive, although by no means all 

49 inclusive. To avoid overlap with the existing excellent reviews [9-11], results reported within the past 

50 few years are featured in words that nonspecialists can readily understand. 

51 2. Background of nanotube buckling research 

52 The term "buckling" means a deformation process in which a structure subjected to high stress un- 

53 dergoes a sudden change in morphology at a critical load [12]. A typical example of buckling may be 

^Together with the three famous nanocarbons, "nanodiamond" deserves attention as the fourth member of nanocarbon 
materials, being classic yet novel. Nanodiamond was initially synthesized by Volkov [5] in 1963, and recently, a wide variety 
of applications has been proposed; see Refs. [6,7] for instance. 

forthcoming monograph [8], which covers broad topics of carbon nanotube deformation, including buckling, will be 
of great help to readers who are interested in the subject. 
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Figure 1. Schematic diagram of buckling of an elastic beam under axial compression: (a) 
pristine beam; (b) axial compression for a small load; (c) buckling observed beyond a critical 
load. 




(a) (b) (c) 



54 observed when pressing opposite edges of a long, thin elastic beam toward one another; see Fig. 1. For 

55 small loads, the beam is compressed in the axial direction while keeping its linear shape [Fig. 1(b)], and 

56 the strain energy is proportional to the square of the axial displacement. Beyond a certain critical load, 

57 however, it suddenly bends archwise [Fig. 1(c)] and the relation between the strain energy and displace- 

58 ments deviates significantly from the square law. Besides axial compression, bending and torsion give 

59 rise to buckling behaviors of elastic beams, where the buckled patterns strongly depend on geometric 

60 and material parameters. More interesting is the elastic buckling of structural pipe-in-pipe cross sections 

61 under hydrostatic pressure [13,14]; pipe-in-pipe (i.e., a pipe inserted inside another pipe) applications 

62 are promising for offshore oil and gas production systems in civil engineering. 

63 The above argument on macroscopic elastic objects encourages to explore what buckled patterns are 

64 obtained in carbon nanotubes. Owing to their nanometric scales, similarities and differences in buck- 

65 led patterns compared with macroscopic counterparts should not be trivial at all. This complexity has 

66 motivated tremendous efforts toward the buckling analysis of carbon nanotubes under diverse loading 

67 conditions: axial compression [15-25], radial compression [26-60], bending [38,61-68], torsion [69- 

68 74], and their certain combinations [75-79]. Such extensive studies have been driven primarily by the 

69 following two facts. One is the excellent geometric reversibility of nanotubes against mechanical de- 

70 formation; that is, their cylindrical shapes are reversible upon unloading without permanent damage to 

71 the atomic structure. In addition, carbon nanotubes exhibit high fatigue resistance; therefore, they are 

72 the promising medium for the mechanical energy storage with extremely high energy density [80]. The 

73 other fact is the substantial influence of buckling on their physical properties. It was recently shown that, 

74 just as one example, carbon nanotubes undergoing an axial buckling instability have potential utility 
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75 as a single-electron transistor [81-83] and can play a crucial role in developing nanoelectromechanical 

76 systems. 

77 Microscopy measurements are powerful means of examining the nonlinear response of nanotubes 

78 against external loading. For instance, atomic force microscopy (AFM) was utilized to reveal the force- 

79 distance curve of nanotubes while buckling [84]. More direct characterizations of nanotube buckling 

80 were obtained by in situ transmission electron microscopy (TEM) [85-87], as partly demonstrated in 

81 Fig. 3 (see §4.2). However, the experimental investigation of nanotube buckling remains a challenge 

82 because of difficulties in manipulation at the nanometric scale. This is a reason why both theoretical 

83 and numerical approaches have played an important role in exploring the buckling behavior of nan- 

84 otubes. In theoretical studies, carbon nanotubes are commonly treated as beams or thin- shell tubes with 

85 certain wall thickness and elastic constants [33,41,88-97]. Such continuum approximations are less 

86 computationally expensive than atomistic approaches; moreover, the obtained formulations can be rela- 

87 tively simple in general. It is noteworthy that, by substituting appropriate values into elastic constants, 

88 continuum-mechanics approaches provide a unified framework [95] that accounts for the critical buck- 

89 ling strains and buckled morphologies under various loading conditions covering compression, bending, 

90 torsion, etc. 

91 3. Resilience and sensitivity to buckling 

92 Special emphasis should be placed on the fact that carbon nanotubes exhibit many intriguing post- 
93 buckling morphologies: Radial corrugations (see §5.2) and surface rippling (§9) are typical examples. 

94 One of the most outstanding features of postbuckled nanotubes is their geometric reversibility upon un- 

95 loading. Indeed, experiments have shown that the buckling deformation can be completely recovered 

96 when the load is removed [61-63,98-101]. The marked structural resilience is primary because of (i) 

97 the large in-plane rigidity of graphene sheets rather than low bending rigidity [102] and (ii) the intrinsic 

98 hollow geometry with extremely large aspect ratio that carbon nanotubes exhibit. It was suggested that 

99 the resilience makes it possible to use the nanotubes as a gas [73,103,104] or water pipeline [105] whose 

100 permeability can be tuned by mechanical deformation. 

101 Apart from the structural resilience, the sensitivity of carbon nanotube properties to buckling is worthy 

102 of attention. In fact, the breakdown of the structural symmetry resulting from the buckling triggers 

103 sudden changes in physical and mechanical properties of nanotubes, including thermal conductivity 

104 reduction [106-108], a radial breathing-mode frequency shift [109], the emergence of interwall sp^ 

105 bondings [110], and electromechanical responses under bending [111-113] and torsion [114,115], to 

106 name a few. In addition, the buckling-induced reduction in nanotube stiffness not only impairs the ability 

107 of nanotubes to sustain external loadings as reinforced fibers in nanocomposites [99,100] but also gives 

108 rise to large uncertainties in the vibration behavior of nanotubes as nanoscale resonators [63,116,117]. 

109 These buckling-property relations can significantly influence the performance of nanotubes as structural 

110 or functional elements, thus implying the need of a huge amount of effort that has been made for the 

1 1 1 study of nanotube buckling. 
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Figure 2. Axially buckled SWNT pattern deduced from molecular dynamics simulations. 
[Left] Upper panel: Energy-strain curve of a (10,10) SWNT with a length of 9.6 nm under 
axial compression. Lower panel: Typical tube geometry (a) before and (b) after buckling, 
respectively. [Right] Upper panel: Curve of a compressed (10,10) SWNT with 29.5 nm 
length. Lower panel: Snapshots of the tube (a) before buckling, (b) after column buckling, 
and (c) undergoing shell buckling. Reprinted from Ref. [118]. 



100 




112 4. Axial compression buckling 

113 4.1. Shell buckling or column buckling ? 

114 Buckled patterns of single- walled nanotubes (SWNTs) under axial compression depend on their as- 

115 pect ratio [15,18,33], which equals the ratio of length to diameter of nanotubes. Roughly, a thick and 

116 short SWNT (i.e., with small aspect ratio) undergoes shell buckling while keeping a straight cylindrical 

117 axis, whereas a thin and long one tends to exhibit a series of shell and column (or Euler) buckling. 

118 The shell buckling process is depicted in the left panel of Fig. 2 [118], where a (10,10) SWNT with 

119 a length of 9.6 nm and an aspect ratio of ^7 was chosen."^ It is seen that the strain energy increases 

120 quadratically with strain at the early prebuckling stage. At a critical strain of 3.5%, a sudden drop in 

121 energy is observed,^ corresponding to the occurrence of shell buckling. During the postbuckling stage, 

122 the strain energy exhibits a linear relationship with strain. The linear growth in energy is understood 

123 by the primary role of the change in carbon-carbon (C-C) bond angles, rather than that of the bond 

124 length variation, in determining the energy-strain relation after buckling. Detailed analyses in Ref. [118] 

^ Recently, shell buckling behavior of short nanotubes with aspect ratior ~ 1 or less was considered, showing the significant 
dependence of the buckling strain on the nanotube length [119]. 

^The critical strain of shell buckling is inversely proportional to the tube diameter [120]. 
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Figure 3. (a)-(f) Series of TEM images of deformation processes for MWNTs initiated by 
applying compressive force in the sample direction, (g) Force-displacement diagram. The 
points indicated by arrows correspond to the TEM images in (d) and (e). Reprinted from 
Ref. [85]. 




125 showed that within the post-buckling regime, the variation in bond angles between neighboring C-C 

126 bonds becomes significant while C-C bond lengths show less variation; this results in an almost constant 

127 axial stress, as deduced from Fig. 2. 

128 With increasing aspect ratio, the buckling mode switches to a column buckling mode owing to the 

129 increased flexibility of the tube. Column buckling is seen in the right panel of Fig. 2 for a much longer 

130 (10,10) SWNT (29.5 nm in length with an aspect ratio of ^22). A sudden drop in strain energy occurs 

131 at a strain level of 1.6%, beyond which the center is displaced in a transverse direction away from its 

132 original cylindrical axis. When the already column-buckled SWNT is further compressed, the structure 

133 curls further and a second drop in strain energy is observed at a strain level of 2.2%. The second drop 

134 corresponds to the onset of another buckling, which is responsible for releasing the excess strain energy. 

135 The tube geometry at this point indicates that the SWNT undergoes shell buckling. With further axial 

136 compression, a linear relationship is observed between energy and strain. This result indicates a column- 

137 to shell-buckling transition of SWNTs with large aspect ratio. ^ 

138 It is important to note that the initial buckling modes, corresponding to the first drop in energy-strain 

139 curve, are different between large- and small-aspect-ratio SWNTs. This fact necessitates an examination 

140 of the validity of continuum-mechanics models for the buckling of SWNTs. Careful assessments of 

141 the continuum approximations have been reported [95-97,120,121], indicating the need to properly use 

^Similar sequential transitions can be observed when a SWNT is bent [65]. See §6.3 for details. 
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142 different models depending on the aspect ratio. As to their consistency with atomistic simulation results, 

143 readers can refer to Ref. [120] in which a list of critical strain data under different conditions is detailed. 



144 4.2. Force-displacement curve 

145 We now turn to experimental facts [122]. Because of the difficulty in sample preparation and ma- 

146 nipulation, only a few attempts have been made to perform axial buckling measurements [85-87]. In 

147 particular, experimental realization of shell buckling under compression has largely behind, though its 

148 signature has been obtained via nanoindentation [123]. Hence in the following discussion, we focus our 

149 attention to the column buckling measurements. 

150 The pioneering work [85] is presented in Fig. 3; The TEM images of (a)-(f) clarify a series of defor- 

151 mation processes for multiwalled nanotubes (MWNTs) initiated by applying a compressive force in the 

152 nearly axial direction.'^ Figure 3(g) shows the corresponding force-displacement diagram. The force at 

153 the initial stage is almost proportional to the displacement [left to the point (d)], indicating the elastic 

154 region, followed by an abrupt decrease at (e). The two points indicated in Fig. 3(g) correspond to the 

155 TEM images in panels (d) and (e), respectively. In Fig. 3(g), the curve right to the point (e) maintains a 

156 slight upward slope. The reason for this post-buckling strength may be due to sequential emergence of 

157 different buckling patterns with increasing the displacement. 

158 A more sweeping measurement on the nanotube resilience was performed for the MWNT with a 

159 higher aspect ratio (^80). Figure 4 shows the resulting force-displacement curve and graphical illus- 

160 tration of the buckling process [84]. An important observation is a negative stiffness region (labeled by 

161 "4" in the plot) that begins abruptly. The sharp drop in force with increasing axial strain, observed at 

162 the boundary of regions (3) and (4), is attributed to the kinking of the MWNT as depicted in the lower 

163 panel. After the kinking takes place, the system is mechanically instable; this behavior is consistent with 

164 the mechanics of kinking described in Refs. [15,85]. The instability seen in region (4) is reproducible 

165 through cyclic compression, which opens up the possibility of harnessing the resilient mechanical prop- 

166 erties of MWNTs for novel composites.^ 

167 5. Radial compression buckling 

168 5.7. Uniaxial collapse ofSWNTs 

169 Radial pressure can yield a distinct class of buckling, reflecting the high flexibility of graphene sheets 

170 in the normal direction. In fact, radial stiffness of an isolated carbon nanotube is much less than ax- 

171 ial stiffness [124], which results in an elastic deformation of the cross section on applying hydrostatic 

172 pressure [26-31,34-37,39,40,42,44] or indentation [56,59,125]. Experimental and theoretical studies, 

173 focused on SWNTs and their bundles, revealed flattening and polygonalization in their cross section 

174 under pressures on the order of a few gigapascals [26,29]. Nevertheless, existing results are rather scat- 

^The actual direction of the applied force deviates from the exact axial one by approximately 15°, as indicated in Fig. 3(a). 
^The instability at (8) and the sharp rise at (9) during unloading stem from the tip pulled out of contact from the sample, 
while the nanotube end remains in contact. 
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Figure 4. [Top] Force-displacement curve of an MWNT with an aspect ratio of ^80 under 
cyclic axial loading. This inset shows a microscopy image. [Bottom] Schematic of the 
change in the MWNT configuration during the buckling process. The labels correspond to 
those indicated in the force-displacement curve. Reprinted from Ref. [84]. 
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175 tered, and we are far away from a unified understanding; for example, the radial stiffnesses of nanotubes 

176 estimated thus far vary by up to three orders of magnitude [31,34,35,40,43,56,59,124]. 

177 The overall scenario of SWNT deformation under hydrostatic pressure is summarized in Fig. 5 [35]. 

178 With increasing pressure, cross sections of SWNTs deform continuously from circular to elliptical, and 

179 finally to peanut-like configurations [35,39,126].^ The radial deformation of carbon nanotubes strongly 

180 affects their physical and structural properties. For instance, it may cause semiconductor-metal tran- 

181 sition [127,128], optical response change [129], and magnetic moment quenching [130] in nanotubes. 

182 From a structural perspective, the radial collapse can give rise to interwall sp^ bonding between adja- 

^For larger radius SWNTs, the peanut-like deformed structure can be transformed to dumbbell-like configurations by van 
der Waals (vdW) attractions between the opposite walls of nanotubes. The latter structure is energetically stable even when 
the applied force is unloaded. 
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Figure 5. Long and short diameters of a (10,10) SWNT as a function of applied hydrostatic 
pressure. The shape of the cross section at some selected pressures is plotted at the bottom 
of the figure. Reprinted from Ref. [35]. 
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183 cent concentric walls [131,132], which may increase nanotube stiffness and therefore be effective for 

184 high-strength reinforced composites [133-136]. 

185 A bundle of nanotubes (i.e., an ensemble of many nanotubes arranged parallel to each other) can 

186 exhibit similar radial collapse patterns to those of an isolated nanotube under hydrostatic pressure. Figure 

187 6 shows [137] the volume change of a bundle of (7,7) SWNTs and a bundle of (12,12) SWNTs as 

188 a function of the applied hydrostatic pressure; the data for a bundle of (7,7)@(12,12) double-walled 

189 nanotubes (DWNTs) is also shown in the same plot. The (12,12) SWNT bundle, for instance, collapses 

190 spontaneously at a critical pressure of 2.4 GPa, across which the cross section transforms into a peanut- 

191 like shape. Two other bundles provide higher critical pressures, as follows from the plot. An interesting 

192 observation is that the transition pressure of the (7,7) tube, which is nearly 7.0 GPa when the tube is 

193 isolated, becomes higher than 10.5 GPa when it is surrounded by the (12,12) tube. This means that 

194 the outer tube acts as a "protection shield" and the inner tube supports the outer one and increases its 

195 structural stability; this interpretation is consistent with the prior optical spectroscopic measurement 

196 [138]. This effect, however, is weakened as the tube radius increases owing to the decreasing radial 

197 stiffness of SWNTs. 



198 5.2. Radial corrugation ofMWNTs 

199 In contrast to the intensive studies on SWNTs (and DWNTs), radial deformation of MWNTs remains 

200 relatively unexplored. Intuitively, the multilayered structure of MWNTs is expected to enhance the radial 

201 stiffness relative to a single-walled counterpart. However, when the number of concentric walls is much 
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Figure 6. Change in the relative volume of the (7,7) @ (12, 12) DWNT bundle and the corre- 
sponding SWNT bundles as a function of hydrostatic pressure. Reprinted from Ref. [137]. 
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202 greater than unity, outside walls have large diameters, so external pressure may lead to a mechanical 

203 instability in the outer walls. This local instability triggers a novel cross-sectional deformation, called 

204 radial corrugation [50], of MWNTs under hydrostatic pressure. 

205 Figures 7(a) and (b) illustrate MWNT cross-sectional views of two typical deformation modes: (a) 

206 elliptic (n = 2) and (b) corrugation (n = 5) modes. In the elliptic mode, all constituent walls are radially 

207 deformed. In contrast, in the corrugation mode, outside walls exhibit significant deformation, whereas 

208 the innermost wall maintains its circular shape. Which mode will be obtained under pressure depends 

209 on the number of walls, A^, and the core tube diameter D of the MWNT considered. In principle, larger 

210 and smaller D favor a corrugation mode with larger n. 

211 Figure 7(c) shows the critical buckling pressure Pc as a function of for various values of D. An 

212 initial increase in Pc at small (except for = 1 .0 nm) is attributed to the enhancement of radial stiffness 

213 of the entire MWNT by encapsulation. This stiffening effect disappears with further increase in A^, 

214 resulting in decay or convergence of Pc{N). A decay in Pc implies that a relatively low pressure becomes 

215 sufficient to produce radial deformation, thus indicating an effective "softening" of the MWNT. The two 

216 contrasting types of behavior, stiffening and softening, are different manifestations of the encapsulation 

217 effect of MWNTs. It is noteworthy that practically synthesized MWNTs often show D larger than those 

218 presented in Fig. 7(c). Hence, Pc{N) of an actual MWNT lies at several hundreds of megapascals, as 

219 estimated from Fig. 7(c). Such a degree of pressure applied to MWNTs is easily accessible in high- 

220 pressure experiments,^^ supporting the feasibility of our theoretical predictions. 

^^A radial pressure large enough to cause corrugation can be achieved by electron-beam irradiation [139]; the self-healing 
nature of eroded carbon walls gives rise to a spontaneous contraction that exerts a high pressure on the inner walls to yield 
their radial corrugation [57]. 
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Figure 7. (a) Cross-sectional views of (a) elliptic (n = 2) and (b) corrugated (n = 5) 
deformation modes. The mode index n indicates the wave number of the deformation mode 
along the circumference, (c) Wall-number dependence of critical pressure Pc. Immediately 
above the original circular cross section of MWNTs gets radially corrugated, (d) Stepwise 
increase in the corrugation mode index n. Reprinted from Refs. [50,58]. 
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221 Figure 7(d) shows the stepwise increases in the corrugation mode index n. For all D, the deformation 

222 mode observed just above Pc abruptly increases from n = 2ton>4ata certain value of A^, followed 

223 by the successive emergence of higher corrugation modes with larger n. These successive transitions in 

224 n at N ^ 1 originate from the mismatch in the radial stiffness of the innermost and outermost walls. 

225 A large discrepancy in the radial stiffness of the inner and outer walls results in a maldistribution of the 

226 deformation amplitudes of concentric walls interacting via vdW forces, which consequently produces an 

227 abrupt change in the observed deformation mode at a certain value of A^. 

228 Other types of radial deformation arise when deviate the interwall spacings of MWNTs from the 

229 vdW equilibrium distance 0.34 nm) [60].^^ The simulations show that the cross sections stabilized at 

^^The authors in Ref. [60] say that this conclusion was motivated by an experimental fact that cross sections of MWNTs 
synthesized in the presence of nitrogen are polygonal shapes rather than circular shapes [140,141]. It was argued that the 
polygonization may result from the interlayer thermal contraction upon cooling or interwall adhesion energy owing to the 
increased interwall commensuration [140]. 
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Figure 8. Cross-sectional views of relaxed MWNTs indexed by (2,8)/(4,16)/. . ./(2n,8n). 
The wall numbers n are 5, 10, 15, 20, and 25 from left to right, and all the MWNTs are 20 
nm long. Reprinted from Ref. [60]. 




230 polygonal or water-drop-like shapes, depending on the artificially expanded interwall spacings.^^ Figure 

231 8 depicts the cross-sectional configurations of relaxed MWNTs. It is seen that the 15-walled tube is 

232 stabilized at a polygonal cross-sectional configuration with six rounded corners. For the 20- and 25- 

233 walled ones, the configuration becomes asymmetric, featuring a water-drop-like morphology. 

234 From an engineering perspective, the tunability of the cross-sectional geometry may be useful in 

235 developing nanotube-based nanofluidic [142-144] or nanoelectrochemical devices [145,146] because 

236 both utilize the hollow cavity within the innermost tube. Another interesting implication is a pressure- 

237 driven change in the quantum transport of tt electrons moving along the radially deformed nanotube. It 

238 has been known that mobile electrons whose motion is confined to a two-dimensional, curved thin layer 

239 behave differently from those on a conventional flat plane because of an effective electromagnetic field 

240 [147-152] that can affect low-energy excitations of the electrons. Associated variations in the electron- 

241 phonon coupling [153] and phononic transport [154] through the deformed nano-carbon materials are 

242 also interesting and relevant to the physics of radially corrugated MWNTs. 

243 6. Bend buckling of SWNTs 

244 6.1. Kink formation 

245 The buckling of SWNTs under bending was pioneered in 1996 [61] using high-resolution electron 

246 microscopes and molecular-dynamics (MD) simulation. Figure 9(a) shows a TEM image of a bent 

247 SWNT with a diameter of 1.2 nm [61]. By bending an initially straight SWNT, its outer and inner sides 

248 undergo stretching and compression, respectively. As a result, it develops a single kink in the middle, 

249 through which the strain energy on the compressed sides is released. Upon removal of the bending 

250 moment, it returns to the initial cylindrical form completely without bond breaking or defects. This 

251 observation clearly proves that SWNTs possess extraordinary structural elastic flexibility. 

^^Prior to structural optimization, the initial cross sections of all the MWNTs are of circular shape and the interwall spacing 
is 0.359 nm, which is 0.19 nm larger than the equilibrium spacing of two graphene sheets (=0.34 nm). 



Version December 22, 201 1 submitted to Materials 



13 of 38 



Figure 9. (a) Kink structure formed in an SWNT with diameters of 1.2 nm under bending. 
The gap between the tip of the kink and the upper wall is about 0.4 nm. (b) Atomic structure 
around the kink reproduced by computer simulations. The shaded circles beneath the tube 
image express the local strain energy at the various atoms, measured relative to a relaxed 
atom in an infinite graphene sheet. The strain energy scale ranges from to 1.2 eV/atom, 
from left to right, (c) Total strain energy (in dimensionless units) of an SWNT of diameter 
^1.2 nm as a function of the bending angle up to 120°. The dip at ^30° in the curve is 
associated with the kink formation. Reprinted from Ref. [61]. 




252 Figure 9(b) presents a computer simulated reproduction of the kink experimentally observed, provid- 

253 ing atomistic and energetic information about the bending process. The overall shape of the kink, along 

254 with the distance of the tip of the kink from the upper wall of the tube, is in quantitative agreement with 

255 the TEM picture of Fig. 9(a). The coding denotes the local strain energy at the various atoms, measured 

256 relative to a relaxed atom in an infinite graphene sheet. In all simulations, the same generic features 

257 appear: Prior to buckling (< 30°), the strain energy increases quadratically with the bending angle, cor- 

258 responding to harmonic deformation [see Fig. 9(c)]. In this harmonic regime, the hexagonal rings on 

259 the tube surface are only slightly strained and the hexagonal carbon network is maintained. Beyond the 

260 critical curvature, the excess strain on the compressed side reaches a maximum and is released through 

261 the formation of a kink, which increases the surface area of the bending side. This is accompanied by a 

262 dip in the energy vs. bending angle curve, as shown in Fig. 9(c). Following the kink formation, the strain 

263 energy increases approximately linearly until bond breaking occurs under quite large deformation. A 

264 similar characteristic energy-strain curve, an initial quadratic curve followed by a linear increase, arises 

265 in the case of axial compression [15], as we learned in §4.2. 

266 6.2. Diameter dependence 

267 The geometrical size is a crucial factor for determining buckling behaviors of SWNTs under bending. 

268 For instance, those with a small diameter can sustain a large bending angle prior to buckling, and vice 

269 versa [155]. Figure 10(a) shows MD simulation data, which show a monotonic increase in the critical 

270 curvature /^c for with a reduction in the nanotube diameter d. The relationship between /^c and d can be 

271 fitted as [15,155] /^c oc rf"^, which holds regardless of the nanotube chirality. 
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Figure 10. (a) Relationship between critical bending buckling curvature and nanotube 
diameter d obtained from MD analyses. The tube length is fixed at 24 nm. (b) The 
length/diameter (L/d) aspect ratio dependence of k^c for the tube chiralities of (5,5), (9,0), 
and (10,10). Reprinted from Ref. [155]. 
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272 In addition to the diameter dependence, the critical curvature of SWNTs is affected by the length/diameter 

273 (L/d) aspect ratio. It follows from Fig. 10(b) that [155] /^c is almost constant for sufficiently long nan- 

274 otubes such that 10 < L/rf < 50 or more, whereas it drops off for short nanotubes satisfying L/d < 10. 

275 The critical aspect ratio that separate the two regions is sensitive to the tube diameter [as implied by the 

276 inset of Fig. 10(b)], but it is almost independent of chirality. 

277 6.3. Transient bending 

278 We have learned in §4 that for relatively thin SWNTs, the buckling is characterized by a disconti- 

279 nuity in the energy curve [Fig. 9(c)]. Across the buckling point, the bending angle dependence of the 

280 deformation energy changes suddenly from being quadratic in the prebuckling regime to linear in the 

281 postbuckling regime. This is, however, not the case for larger diameter nanotubes. As the diameter is 

282 increased, a second discontinuity appears in the strain-energy curve at a larger bending angle than the 

283 first one [65]. The origin of the two discontinuities can be accounted for by inspection of Fig. 11. The 

284 bottom image [Fig. 1 1(b)] shows a just-buckled wall of a (30, 30) SWNT corresponding to the first dis- 

285 continuity in the energy curve, in which the buckled side is far from the opposite side. Therefore, more 

286 bending is required to bring the two sides close enough, as observed in Fig. 11(a), which results in the 

287 second discontinuity. 

288 The thick-nanotube's buckling behavior mentioned above is illustrated in Fig. 12(a), where the defor- 

289 mation energy U is plotted as a function of bending angle 9 for a (30, 30) SWNT [65]. Three distinct 

290 deformation regimes are observed, clearly separated by two discontinuities at 6^ = 12° and 32°. In the 
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Figure 11. Predicted shape of SWNTs just after buckling, based on MD simulations for (a) 
a 15.7-nm-long (10, 10) SWNT at the bending angle 9 = 43° and (b) a 23.6-nm-long (30, 
30) SWNT at e = 23°. Note the difference in scale. Reprinted from Ref. [65]. 




291 initial elastic regime, U exhibits a quadratic dependence on 9, whereas the cross section experiences 

292 progressive ovalization as the bending angle increases, culminating to the shape in Fig. 12(b). The buck- 

293 ling event is marked by an abrupt transition from the oval cross section to one with the flat top shown in 

294 Fig. 12(c). As the bending angle increases beyond the first discontinuity [i.e., during the transient bend- 

295 ing regime (TBR) indicated in Fig. 12(a)], the flat portion of the top wall expands continuously across 

296 the nanotube [Figs. 12(d) and (e)].^^ As a result, the top-to-bottom wall distance decreases gradually, re- 

297 ducing the tube cross section at the buckling site. The deformation energy curve in the TBR is no longer 

298 quadratic; in fact, the exponent becomes less than unity. When the approaching opposite walls reach the 

299 vdW equilibrium distance of 0.34 nm [Fig. 12(f)], the cross section collapses, forming the kink, and a 

300 second discontinuity is observed. 

301 7. Bend buckling of MWNTs 

302 7.1. Emergence of ripples 

303 Following the discussions on SWNTs in §6, we look into the bend buckling of MWNTs. The differ- 

304 ence in the mechanical responses between SWNTs and MWNTs lies in the presence of vdW interactions 

305 between the constituent carbon layers. Apparently, thicker MWNTs with tens of concentric walls seem 

306 stiff er than few- walled, thin MWNTs against bending, since the inner tubes of MWNTs may reinforce 

307 the outer tubes via the vdW interaction. However, the contrary occurs. In fact, whereas MWNTs with 

308 small diameter exhibited a bending stiffness of around 1 TPa, those with larger diameter were much 

309 more compliant, with a stiffness of around 0.1 TPa [63]. This dramatic reduction in the bending stiffness 

310 was attributed to the so-called rippling effect, i.e., the emergence of a wavelike distortion on the inner 

311 arc of the bent nanotube [99,100,156-166]. 

^^Interestingly, a buckled SWNT in the TBR is fully reversible. If bending is stopped before the second discontinuity 
occurs, unbending recovers the cross-sectional shapes at the buckling point. 
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Figure 12. [Left] Deformation energy U for a 23.6-nm-long (30,30) SWNT as a function 
of the bending angle 9. The symbols a-e attached to the curve indicate the points for which 
the tube shape and cross section at the buckling point are shown in the images (b)-(f) on the 
right. "TBR" denotes the transient bending regime. Reprinted from Ref. [65]. 




312 Figure 13(a) presents a clear example of the rippled MWNT structure [99]. The tube diameter is ^31 

313 nm and it is subjected to a radius of bending curvature of ^400 nm. Enhanced images at the ripple region 

314 [63] are also displayed in Figs. 13(b)-(d), though they are not identical to the specimen in Fig. 13(a). 

315 The amplitude of the ripple increased gradually from inner to outer walls, being essentially zero for the 

316 innermost core tube to about 2 to 3 nm for the outermost wall. Such rippling deformation induces a 

317 significant reduction in the bending modulus, as has been explained theoretically by solving nonlinear 

31 8 differential equations [ 1 67] . 

319 7.2. Yoshimura pattern 

320 Precise information about the membrane profile and the energetics of the rippling deformation, which 

321 are unavailable in experiments, can be extracted from large-scale computer simulations. Figure 14(a) 

322 shows a longitudinal cross section of the equilibrium configuration [161]. This image is the compu- 

323 tational analog of the TEM slices of rippled, thick nanotubes reported in the literature [63,160]. The 

324 simulations reproduce very well the general features of the observed rippled nanotubes: nearly periodic 

325 wavelike distortions, whose amplitudes vanish for the inner tubes and smoothly increase toward the outer 

326 layer. 
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Figure 13. (a) Under high bending, MWNTs form kinks on the internal (compression) 
side of the bend, (b)-(d) High-resolution TEM image of a bent nanotube (with a radius of 
curvature of ^400 nm), showing the characteristic wavelike distortion. The amplitude of the 
ripples increases continuously from the center of the tube to the outer layers of the inner arc 
of the bend. Reprinted from Refs. [63,99]. 



(a) 




327 A remarkable finding in the simulations is that the rippling deformation closely resembles the Yoshimura 

328 pattern^"^ (a diamond buckling pattern). We can see that the rippling profile in Fig. 14(b) consists not of 

329 a simple linear sequence of kinks but of a diamond-like configuration of kinks on the compressed side. 

330 Such a Yoshimura pattern is well known as characterizing the postbuckling behavior of cylindrical elas- 

331 tic shells on a conventional macroscopic scale. The pattern has the interesting geometric property of 

332 being a nearly isometric mapping of the undeformed surface, at the expense of creating sharp ridges and 

333 furrows. 

334 The rippling deformation, peculiar to thick MWNTs, is a consequence of the interplay between the 

335 strain-energy relaxation and the vdW energy increment. As intuitively understood, the low bending 

336 rigidity of individual graphitic sheets, relative to their large in-plane stiffness, makes it possible to release 

337 effectively a significant amount of the membrane strain energy at the expense of slight flexural energy. 

338 As a result, rippled MWNTs have a significantly lower strain energy than uniformly bent MWNTs. 

339 Figure 15 shows the energy of bent MWNTs as a function of the bending curvature for a 34- walled 

340 nanotube [161]. When the nanotube is uniformly bent, the strain energy grows quadratically with re- 

341 spect to the curvature. For such a uniform bending, the vdW energy gives almost no contribution to 

342 deformation, and therefore, the total energy also follows a quadratic law. However, the actual behavior 

343 of the system greatly deviates from this ideal linearly elastic response. As can be observed in Fig. 15, 

^"^The Yoshimura pattern is a special kind of surface deformation mode occurring in thin- walled cylindrical shells subjected 
to large lateral load. It is named after Prof. Yoshimura [168], a Japanese theoretician of fracture mechanics, and its profile is 
characterized by a periodic diamond-like corrugation. 
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Figure 14. Rippling of a 34-walled carbon nanotube: (a) longitudinal section of the central 
part of the simulated nanotube and (b) the morphology of the rippled MWNT reminiscent 
of the Yoshimura pattern. Highlighted are the ridges and furrows, as well as the trace of the 
longitudinal section. Reprinted from Ref. [161]. 
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344 the rippling deformation leads to much lower values of strain energy and an increase in vdW energy. 

345 The evolution of the total energy E'tot with respect to curvature radius R is very accurately fitted by 

346 E'tot oc with a = 1.66; this response differs from that predicted by atomistic simulations of SWNTs 

347 or small, hollow MWNTs, both of which exhibit an initial quadratic growth (a = 2) in the elastic regime, 

348 followed by a linear growth (a = 1) in the postbuckling regime. The results in Figs. 14 and 15 evidence 

349 the failure of the linear elasticity and linearized stability analysis to explain the observed well-defined 

350 postbuckling behavior (1 < a < 2) for thick nanotubes, implying the need for a new theoretical frame- 

351 work based on nonlinear mechanics. 

352 8. Twist buckling 

353 8.1. Asymmetric response of SWNTs 

354 Similarly to bending situations, SWNTs under torsion exhibit a sudden morphological change at a 

355 critical torque, transforming into a straight-axis helical shape. The crucial difference from bending 

356 cases is that, under torsion, the critical buckling torque of SWNTs depends on the loading direction, 

357 i.e., whether the tube is twisted in a right-handed or left-handed manner [171,172]. This load-direction 

358 dependence originates from the tube chirality, which breaks the rotational symmetry about the tube axis. 

359 For example, the twisting failure strain of chiral SWNTs in one rotational direction may even be 25% 

360 lower than that in the opposite direction [172]. Moreover, symmetry breaking causes coupling between 

361 axial tension and torsion, giving rise to an axial- strain-induced torsion of chiral SWNTs [173]. This 

^^Interestingly, the energetics of MWNT bucking under bending are altered by inserting cross-linking (i.e., sp^ covalent 
bonding) between adjacent walls, as presented in Refs. [169,170]. 
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Figure 15. Energy curves for a bent 34-walled nanotube with respect to the bending curva- 
ture. Shown are the strain energy for fictitiously uniform bending (squares), the strain energy 
for actually rippled deformation (crosses), and the total energy (i.e., sum of the strain energy 
and the vdW one) for rippled deformation (circles). Reprinted from Ref. [161]. 



0.05 




362 intriguing coupling effect shows that a chiral S WNT can convert motion between rotation and translation, 

363 thus promising a potential utility of chiral SWNTs as electromechanical device components. 

364 The effect of structural details on buckling of a torsional SWNT was explored using MD calculations 

365 [171]. Figure 16 shows morphology changes of an (8,3) SWNT under torsion. Its torsional deformation 

366 depends significantly on the loading direction. Under right-handed rotation, the tube buckles at a critical 

367 buckling strain of ^7.6%, which is significantly larger than that (^4.3%) under left-handed rotation. 

368 Figure 17 summarizes a systematic computation [171] of the critical buckling strains in both twisting 

369 (7cr) and untwisting (7ci) directions as a function of tube chiral angle; a loading-direction-dependent 

370 torsional response of chiral tubes is clearly observed. Special attention should be paid to the fact that the 

371 maximum difference between and %\ is up to 85%. This clear difference in the mechanical response 

372 suggests particular caution in the use of carbon nanotubes as torsional components (e.g., oscillators and 

373 springs) of nanomechanical devices [1 14,1 17,174,175]. 

374 8.2. Nontrivial response ofMWNTs 

375 In contrast to SWNTs cases, fewer studies on MWNTs under torsion were reported because of their 

376 complex structures and computational costs. For twisted DWNTs, MD simulations have revealed [72,74] 

377 a nontrivial buckling mode involving a few thin, local rims on the outer tube while the inner tube shows 

378 a helically aligned buckling mode [Figs. 18(a)-(c)]. These distinct buckling modes of the two concentric 

379 tubes imply that a conventional continuum approximation in which it is postulated that the buckling 
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Figure 16. Morphological changes for a (8, 3) nanotube under torsion. Applied strain and 
its direction are indicated beneath the diagram; the digit 0.05, for example, corresponds 
to the strain of 5%, and the sign + (— ) indicates right(left)-handed rotation. Under right- 
handed rotation, the tube buckles at a critical buckling strain =7.6%, whereas it buckles 
at 7ci = 4.3% under left-handed rotation. Reprinted from Ref. [171]. 
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380 modes of all the constituent tubes have the same shape fails for analyzing the torsional responses of 

381 DWNTs. 

382 When increasing the number of constituent walls to far more than two, we acquire torsional rippling 

383 deformations [176].^^ It was numerically found that [165] the amplitude in the torsional rippling of 

384 MWNTs can be accurately described by a simple sinusoidal shape function, as confirmed by Figs. 18(d) 

385 and (f). It is noteworthy that the characteristics of the helical rippling morphology in twisted MWNTs are 

386 different from those in bent MWNTs, i.e., the so-called Yoshimura or diamond buckling pattern. Struc- 

387 turally, torsion-induced rippling is distributed more periodically and uniformly along the tube whereas 

388 bending-induced rippling is located only in the compressive side [161,163]. Energetically, in twisted 

389 MWNTs, high strain energy is stored along the ridge regions, whereas in bent MWNTs, the strain en- 

390 ergy is equally concentrated both at the ridges and furrows. 

391 9. Universal scaling laws under bending and torsion 

392 Interestingly, buckled MWNTs under bending and torsion were found to obey universal scaling laws 

393 that consist of two distinct power-law regimes in the energy-deflection relation [163]. Figure 19 shows 

394 the mechanical response of MWNTs under bending and torsion; the strain energy E vs. bending curva- 

395 ture tv or twisting angle 6 is plotted with the increasing stepwise number of walls from 10 to 40. All the 

396 tested MWNTs exhibit two distinct power-law regimes: a harmonic deformation regime characterized 

^^During torsional rippling, the innermost core tube stores a very high strain energy despite the nearly zero rippling am- 
plitude. Owing to the strong confinement, the strain energy in this layer cannot be released via rippling, which may lead to 
bond breaking and subsequent brittle cracking when the torsional deformation continues to increase [162,164]. 
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Figure 17. Critical buckling shear strains as a function of tube chirality. Some additional 
data for SWNTs with slightly larger or smaller diameters are also presented for reference. 
Reprinted from Ref. [171]. 
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397 by the exponent a = 2 for the relation oc a^:^ or oc (indicated by blue lines in Fig. 19) and an 

398 anharmonic, postbuckling regime with exponent a ^ 1.4 for bending and a ^ 1.6 for torsion (red lines). 

399 The latter nonlinear response corresponds to the ripples of the graphene walls discussed earlier. 

Figures 19(b) and (e) show the data collapse for all the tested nanotubes upon a universal scaling 
law. This law is described by the anharmonic exponent a and the characteristic length scale ^cr- In twist 
situations, it is defined by ^cr = Qc^R^ /L, where L and R are the length and outer radius, respectively 
of the MWNTs considered, and ©cr denotes the critical twisting angle at which the buckling arises. In 
bending cases, ^cr = t^crR^ with /^cr being the critical buckling curvature. Then, the unified law plotted 
in red and blue in Figs. 19(b) and (e) is 

E{x) f {^Rf for l^^l < 4r, 

—r- =^ { (1) 

^ [ for \xR\ > 4r, 

400 with X = QR/L or X = tvR. The actual value of 4r was evaluated as 4r ^0.1 nm for both bending 

401 and twisting cases. It should be emphasized that since 4r has dimension of length, the unified law is size 

402 dependent; for instance, the thicker the MWNTs, the smaller will be the obtained Bcr or /^cr- 

403 10. Challenge and future directions 

404 In this article, I have provided a bird's-eye view on the current state of knowledge on the buckling 

405 properties of carbon nanotubes. The understanding remains far from complete, but new experiments and 

406 theoretical work will no doubt give us a more complete picture and exciting times for both basic and 
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Figure 18. (a)-(c) Helical buckling of (5,0)@(14,0) DWNTs withlengths (a) L = 1.095 nm, 
(b) 4.45 nm, and (c) 6.97 nm. For each tube length, the inner wall shows ordinary buckled 
patterns, whereas the outer wall exhibits nontrivial buckling modes associated with local 
rims, (d)-(f) Helical rippling deformation of a 10-walled nanotube (5,5) @. . . @ (50,50) of 
34 nm in length and 3.4 nm in radius, (d) Longitudinal view, (e) Cross-sectional view, (f) 
Deformation map, with green for ridges and blue for furrows. Reprinted from Refs. [72,165]. 




407 applied research in the realm of nanoscale. Described below are only a few examples of challenging 

408 subjects that may trigger innovation in the nanotube research community. 

409 70.7. Buckling effects on heat transport 

410 Carbon nanotubes demonstrate the excellent thermal conductivity among any known material. When 

411 a carbon nanotube is buckled, however, the localized structural deformation can prohibit ballistic heat 

41 2 transport along the nanotube axis [108]. This results in the decreasing behavior of thermal conductivity of 

413 nanotubes under compressive stress, which is attributed to the increase in the phonon-phonon scattering 

414 rate. Such the buckling-induced reduction in thermal conductivity has important implications of heat 

415 management of nano-scale electronic devices, including the dynamic control of thermal transport and 

416 energy conversion. In view of materials sciences, it is also important to make clear the effect of buckling 

417 on the thermal properties of carbon nanotube composites [177]. Despite of great potential, very limited 

418 number of studies has been conducted on the issue thus far. 

419 From an academic standpoint, the buckling-induced change in the thermal transport poses questions 

420 on the feasibility of conventional heat conduction theory for macroscopic solids. It has been clarified that 

421 the low-dimensional nature of carbon nanotubes gives rise to various intriguing phenomena: the well- 

422 known examples are the lattice soliton based energy transfer [178] and robust heat transport in deformed 

423 nanotubes [179]. These phenomena are beyond the classical way of understanding; therefore, it is inter- 
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Figure 19. (a,d)* Strain energy curves as a function of the bending curvature and the 
twisting angle B. (b,e): Data collapse upon appropriate rescaling. The power-law fits with 
exponents 2 (blue) and a(< 2) (red) are shown for illustration. In all four plots, the number 
of walls increases stepwise from 10 (circles) to 40 (crosses), (c): 40- walled nanotube in pure 
bending, (f): 35-walled nanotube in torsion. The latter two panels present the deformed 
shape (top), Gaussian curvature map (middle, with green being zero, red being positive, and 
blue being negative), and energy density map (bottom, with red being high and blue being 
low). Reprinted from Ref. [163]. 
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424 esting to consider how each class of nanotube buckling (compressive, bending, torsional, etc.) affects 

425 the nontrivial heat transport observed in nanotubes. The results obtained will shed light on unexplored 

426 problems of thermal conduction in carbon nanotubes and related materials. 

427 70.2. Role of defects and imperfections 



428 Carbon nanotubes obtained by practical synthesis possess various kinds of defects [8] such as miss- 

429 ing atoms (called "vacancy defects"), carbon rings other than usual hexagonal ones ("Stone-Wales de- 

430 fects"), or sp^ bonds instead of usual sp^ bonds ("re-hybridization defects"). Hence, understanding the 

431 effect of defects on the mechanical properties of nanotubes is essential in the design of nanotube-based 

432 applications. It has been found that these defects can affect considerably the mechanical strength and 

433 post-elongated morphology of carbon nanotubes. Comprehensive studies on the influence of defects on 

434 their buckling behavior remained, however, lacking in the literature until recently [11]. 
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435 MD simulations performed in the past few years [180-183] suggested that the presence of defects, 

436 particularly one- or two-atom vacancies, may cause a significant reduction in the buckling capacity. An 

437 interesting observation is that the degree of reduction is strongly dependent on the chirality and tem- 

438 perature. For example, in torsional buckling at low temperature, armchair nanotubes are less sensitive 

439 to the presence of defects when compared with their zigzag counterparts [182]. Still theoretical in- 

440 vestigations have been limited to MD simulations; more accurate and quantitative research based on 

441 density-functional methods, for instance, would be desired. Experiments on the defect-induced variance 

442 in the buckling behavior are of course to be addressed in future. 

443 Another interesting subject is to employ the presence of defects as regulator of heat conduction 

444 through carbon nanotubes. With the increase of number of defects, the thermal conductance of nanotubes 

445 rapidly decreases. The reason for this large reduction is that high-frequency phonons which contribute 

446 to thermal transport is strongly scattered by the structural defects. It has been numerically predicted that 

447 [184] even a few structural defects in nanotubes can lead to a strong suppression of thermal transport by 

448 one order of magnitude. This result implies that the structural defects can offer an effective method of 

449 tuning thermal transport of carbon nanotubes; such the tuning of heat transport is advantageous in the 

450 sense that lattice defects can be well controlled during the growth or by irradiations. 

451 10.3. Relevance to chemical reaction 

452 Applications of carbon nanotubes in composite materials often require an understanding and control 

453 of the chemistry and chemical reactivity of the nanotubes' sidewall. This is because the carbon-carbon 

454 bonding state on the outermost graphitic surfaces determines where the chemically sensitive reactions 

455 take place and how the reactions affect the physical properties of carbon nanotubes. For instance, the 

456 functionalization and/or chemisorption on the sidewall of nanotubes enable to increase linking between 

457 nanotubes as filler and a surrounding matrix. Besides, reactivity control of the nanotube sidewall leads 

458 to a novel technique for chemical sensors and drug delivery systems. 

459 An important finding in the context of buckling is that the chemical reactivity of carbon nanotube is 

460 dependent on local surface curvature of the outermost sidewall. It has been theoretically demonstrated 

461 that [185,186] the reactivity of a nanotube is governed by the local atomic structure of carbon atoms on 

462 which the chemisorbing species or functional group can react and/or a stable bond. This results imply 

463 the control of sidewall reactivity by artificial deformation; that is, local chemical reactivity of carbon 

464 nanotubes can be promoted locally by inducing mechanical deformation like buckling. 

465 It should be reminded that the basic concept of the curvature-dependent reactivity was initially pro- 

466 posed in 1993 [187]; almost twenty years have passed since then. Nevertheless, research progress on 

467 the subject seems rather behind [188,189]. Especially, experimental evidence of the proposed reactivity 

468 promotion has been scarce, which may be due to difficulty in precise manipulation of nanotube buckling 

469 and accurate measurement of adsorption capability on the nanotube sidewall. Considering the impor- 

470 tance in view of material science, breakthrough in the chemistry of the deformation-driven reactivity and 

471 its possible application is strongly expected. 
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Figure 20. TEM image of nanotubes and sketches of each nanotube's cross sections for (a) 
a five-wall nanotube with a diameter of 6.7 nm, (b) a double-wall nanotube with a diameter 
of 5.5 nm, and (c) a seven- wall nanotube with a diameter of 6.5 nm, which has the smallest 
hollow diameter (^2.2 nm) among the three specimens. Electron diffraction patterns show- 
ing (d) the superposition of three sets of {hkO} spots taken from a seven-wall nanotube and 
(e) the superposition of four sets of {hkO} spots from a nine-wall nanotube. Reprinted from 
Ref. [3]. 




472 Appendix: 

473 A. Who discovered carbon nanotubes first? 

474 In Appendix, we take a look at the history of who discovered nanotubes and when,^^ for mainly 

475 pedagogical reasons. The answer seems obvious at a glance, but it is actually quite profound. The author 

476 believes that the argument should be of help for younger readers who may be learning the facts for the 

477 first time. 

478 A. 7. lijima's nanotube in 1991 

479 As is well known, most academic and popular literature attributes the discovery of carbon nanotubes 

480 to Sumio lijima of NEC in 1991. Without doubt, it is lijima's seminal article in 1991 [3], entitled 

481 "Helical microtubules of graphitic carbon," that triggered the explosion of carbon nanotube research still 

^^The relevant literature has been probed with great care by Boehm in 1997 [190] and Monthioux and Kuznetsov in 2006 
[191]. The first half of this Appendix is based on the two excellent reviews. 
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Figure 21. The earliest TEM images of carbon nanotubes published in Ref. [192]. 




482 fascinating us today. lijima studied the arc evaporation process that efficiently produced fuUerene (Ceo) 

483 molecules. When analyzing a by-product in the arc evaporation, lijima found large amounts of MWNTs 

484 mixed with faceted graphitic particles; see Fig. 20 for the image obtained. 

485 From a historical perspective, however, lijima's finding is not the first reported carbon nanotube. 

486 Careful analysis of the literature shows that there had been many precedents prior to 1991, in which the 

487 presence of analogous (and almost identical) nanostructures including "carbon tubes" [192] and "hollow 

488 carbon fibers" [193] was uncovered. 

489 A. 2. Carbon nanotubes prior to 1991 

490 The first evidence for nanosized carbon tubes is believed to have been published in 1952 in a Soviet 

491 academic journal [192], almost 40 years before lijima's paper. The article, written by Radushkevich and 

492 Lukyanovich, demonstrated clear TEM images of tubes (diameter ^ 50 nm) made of carbon. Figure 21 is 

493 a reprint of the image given in Ref. [192], which clearly shows carbon filaments exhibiting a continuous 

494 inner cavity. However, this discovery was largely unnoticed as the article was published in the Russian 

495 language. 

496 Subsequently, many other reports followed this Soviet work. For instance, Oberlin, Endo, and 

497 Koyama fabricated in 1976 hollow carbon fibers with nanometer-scale diameters using a vapor- growth 

498 technique [193]. In 1979, Abrahamson presented evidence of carbon nanotubes at a conference paper 

499 together with their characterization and hypotheses for a growth mechanism [194]. In 1981, a group 

500 of Soviet scientists produced carbon tubular crystals and identified them with graphene layers rolled up 

501 into cylinders. Moreover, they speculated that by rolling graphene layers into a cylinder, many different 

502 arrangements of graphene hexagonal nets, such as an armchair-type and a chiral one, are possible. 

Despite the omission of a scale bar, the magnification value indicated in the article allows one to estimate that the 
diameters of the carbon tubes imaged are in the range of 50 nm, i.e., definitely nanosized. 
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503 A. 3. Closing remarks 

504 In summary, the credit for discovering carbon nanotubes should go to Radushkevich and Lukyanovich. 

505 It must be emphasized, however, that merely looking at an object with a microscope is completely dif- 

506 ferent from identifying its structure. That is, lijima not only observed the nanotubes with an electron 

507 microscope but also accurately elucidated its structure from electron diffraction images, and his work 

508 was thus on a different level to preceding studies. 

509 Let us return to Figs. 20(a)-(c). The photographs shown here are the first electron micrographs of 

510 carbon nanotubes captured by lijima. Surprisingly, a multiwalled structure with an interlay er spacing 

511 of less than 1 nm is clearly observed in this photograph. Even more surprisingly, he measured the 

512 electron diffraction shape [Figs. 20(d) and (e)] from a single carbon nanotube. The diffraction strength 

513 from a crystal of a light element (carbon) with a thickness of only several nanometers is extremely 

514 low, and so, this experiment is not easy to perform even nowadays using the best electron microscopes 

515 available. lijima must have had first-class microscope operating skills to make these measurements in 

516 1991, although he did it with ease. It is also astounding that lijima identified this substance, which only 

517 appeared as a thin, needle-like crystal, as a "tube-shaped substance with a helical structure" based on its 

518 electron diffraction shape. How many people in the world would have been able to conclude that it is a 

519 nanotube with a helical structure by looking at the electron diffraction shape? Now that we already know 

520 the structure of carbon nanotubes, it is not so difficult to read the electron diffraction shape. However, 

521 at that time, much experience and sharp insight must have been required by the person observing this 

522 electron diffraction shape for the first time in the world to interpret it correctly as a tube-shaped substance 

523 with a helical structure. 

524 So, who actually discovered carbon nanotubes? To answer this question, first we need to define the 

525 meaning of "discovery of nanotubes." It is probably correct to state that Radushkevich and Lukyanovich 

526 were the first to take photographs of the overall picture of the cylindrical tube- shaped nanocarbon sub- 

527 stance (although photographs taken by other researchers could still be discovered in future). However, 

528 lijima was the first to reveal that the microscopic structure unique to nanotubes was multiwalled and 

529 helical, by using "the eyes" of electron diffraction. The rest is up to the philosophy and preference of 

530 the person who is judging. Whatever the judgment, we should never forget that the knowledge and ben- 

531 efits that we now enjoy in nanotube research are the result of the hard work and achievements of our 

532 predecessors. 
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